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Abstract 

An effective action technique for the time evolution of a closed system consisting 
of one or more mean fields interacting with their quantum fluctuations is presented. 
By marrying large N expansion methods to the Schwinger-Keldysh closed time path 
(CTP) formulation of the quantum effective action, causality of the resulting equa- 
tions of motion is ensured and a systematic, energy conserving and gauge invariant 
expansion about the quasi-classical mean field(s) in powers of 1/N developed. The 
general method is exposed in two specific examples, 0(N) symmetric scalar A$ 4 the- 
ory and Quantum Electrodynamics (QED) with N fermion fields. The A$ 4 case is 
well suited to the numerical study of the real time dynamics of phase transitions 
characterized by a scalar order parameter. In QED the technique may be used to 
study the quantum non-equilibrium effects of pair creation in strong electric fields 
and the scattering and transport processes in a relativistic e + e~ plasma. A simple 
renormalization scheme that makes practical the numerical solution of the equations 
of motion of these and other field theories is described. 



1 Introduction 



The derivation of macroscopic dissipative behavior from fundamental time reversal 
invariant dynamics is a subject at least as old as Boltzmann. In the classical approach 
to this issue two ingredients are necessary. First, non-equilibrium initial conditions for 
the system are considered, and second, some kind of coarse-graining or averaging pro- 
cedure over unobserved degrees of freedom is introduced [1]. The classical prototype 
for this analysis is of course, Boltzmann's equation for the evolution of the molecular 
distribution function in a gas due to collisions between the molecules. Another classic 
example is the Brownian motion of a heavy particle in a fluctuating medium, where 
a clean separation between the "system" and "environmental" degrees of freedom 
to be averaged is assumed in the first step. In both of these examples from clas- 
sical statistical mechanics the phenomenological point of view is the dominant one, 
and there is no question about how entropy growth or loss of information arises. It 
arises either by explicitly replacing the underlying time reversal invariant dynamics 
by stochastic assumptions, or by the separation into "system" and "environment" 
with all detailed information about the latter lost in the averaging procedure (save 
for a few parameters, such as temperature). 

In recent years there has been considerable growth in interest in the study of 
dissipative effects in specifically quantum mechanical systems. Driven largely by the 
impressive progress in the fabrication and control of sensitive microdevices, such as 
tunnel junctions, issues of quantum dissipation have been studied extensively by a 
variety of methods [2]. At the other end of the distance scale, the subject of loss of 
coherence in closed quantum systems has become of interest in quantum cosmology, 
where one would like to understand precisely how the quasi-classical universe we 
observe emerged from (presumably) purely quantal initial conditions [3]. Because of 
their technical simplicity and applicability to these widely different situations, model 
theories of quantum Brownian motion have often served as prototypes in studies of 
dissipation and decoherence. In fact, almost all existing treatments of non-equilibrium 
dissipative phenomena in quantum systems focus on some particular model or class of 
models, such as the coupled oscillators of Caldeira and Leggett [4], chosen primarily 
for the purpose of illustrating some definite feature of the analysis in the simplest way, 
or for describing a particular existing phenomenology. Since relativistic effects are 
unimportant in condensed matter applications, they have been generally neglected in 
those phenomenological models. 

There is no essential problem in extending the analyses of non-equilibrium and 
irreversible processes to systems with both quantum and relativistic features, i.e., to 
relativistic quantum field theories. Such an analysis is not necessary for the primary 
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application of field theories to experimental situations, namely scattering experiments 
involving a few particles, and this explains why it has not been much discussed in 
the literature. However, implicitly assumed in the calculation of a scattering matrix 
or cross section is that the particles in both the distant past (the in state) and 
in the remote future (the out state) of the scattering event are to be treated as 
essentially free particles. It is only this quite restrictive assumption that justifies 
the usual Feynman boundary conditions on Green's functions, and the corresponding 
path integral representation of the scattering matrix between in and out states. 

As successful as this scattering theory has been, it is clear that it is capable of 
addressing only a very limited subclass of questions one might hope to address in 
quantum field theory. The coupling may not be weak enough to justify a perturba- 
tive expansion, there may be no approximately free asymptotic scattering states, or 
we may be interested in more than just scattering probabilities, and of course, such a 
treatment based on asymptotic states cannot teach us anything directly about dissi- 
pation or loss of quantum coherence. Only relatively recently has it become apparent 
that there are many situations of physical interest requiring a detailed description of 
the time- dependent dynamical evolution, where the perturbative scattering formal- 
ism of quantum field theory is wholly inadequate. For example, this is the case for 
strong field electrodynamics in astrophysical plasmas, in nuclear collisions of heavy 
ions where it is possible that a phase transition to the quark-gluon plasma may be 
produced, in phase transitions in the hot, dense early universe, and in the dynamics of 
phase transitions generally. All of these problems require a detailed knowledge of the 
dynamical time evolution of the field configurations, as well as their non-equilibrium 
transport and energy-momentum flow characteristics. It is at that point also that the 
fundamental issue of dissipation in the time reversal invariant physics described by 
the underlying quantum fields asserts itself. 

Despite the apparent need for a time-dependent (and therefore not explicitly co- 
variant) formulation of quantum field theory, to address both the fundamental issues 
and the many interesting applications, the technical complications of the renormaliza- 
tion procedure in a non-covariant treatment have usually been sufficient to dissuade 
all but the most persistent from the task. When renormalization issues are faced 
squarely, it is usually by rather formal methods, such as dimensional continuation. 
Although elegant theoretically this formal regularization technique is not adaptable 
to a computer, to which it should be clear from the outset one will have to resort to 
in order to solve the coupled non-linear partial differential equations of any realistic 
theory. In fact, it is only the advent of modern supercomputers that makes it even 
conceivable to carry out the program outlined in this paper and obtain useful results 
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in realistic physical systems. 

Motivated by these various considerations we wish in this paper to provide a 
fully quantum field theoretic framework for the systematic inclusion of the effects of 
fluctuations on the time evolution of a closed system in a formulation practical for nu- 
merical methods. In terms of our original discussion of "system" and "environment," 
the mean field (or fields) will play the role of the macroscopic subsystem, while the 
fluctuating degrees of freedom become the environment or bath in which the system 
moves. The important difference in our approach will be that the latter will not be 
introduced externally, but rather determined from the correlation functions of the 
theory self-consistently in terms of the time evolving mean field(s). In other words 
we wish to discuss the quantum evolution of a closed system, where no averaging 
over environmental degrees of freedom is ever performed. It is clear that the exact 
solution of this problem could never allow for any dissipation or decoherence, since 
well-defined Hamiltonian and unitary time evolution operators exist in the underly- 
ing field theory. However, we cannot hope to solve the Schwinger-Dyson equations 
of a non-trivial interacting quantum field theory exactly, so that some systematic 
approximation scheme is certainly required. In statistical mechanics the natural ap- 
proximation scheme is to truncate the infinite hierarchy of correlation functions at 
some finite order. It is at this point when higher order correlations are neglected 
that information about the exact unitary evolution may be lost, and effective time 
irreversibilty may enter. This interesting point has been emphasized by Calzetta and 
Hu [5]. 

The problem of the time evolution of a closed system together with its self- 
generated quantum fluctuations is clearly very general, and arises in a number of 
contexts in very different areas of physics. Original motivations in microscopic and 
mesoscopic quantum devices have been mentioned [2], but there are many other pos- 
sible applications. To consider but a few examples, in electromagnetic plasmas the 
classical phenomenon of radiation reaction has its quantum counterpart when particle 
creation and virtual processes are considered. These lead to the damping of plasmon 
modes in strong fields, such as those in the vicinity of rapidly rotating neutron stars 
[6]. In heavy ion collisions of ultrarelativistic nuclei a strongly interacting quark-gluon 
plasma phase may be produced [7]. The theoretical understanding of nonequilibrium 
effects in such a plasma is very rudimentary at present, but is clearly required to 
interpret the data soon to be generated at accelerators such as RHIC or the LHC. 
Similar effects are expected in the real time dynamics of any phase transition, whether 
occuring in terrestrial condensed matter systems or in the early universe. In black 
hole evaporation, the effect of fluctuating degrees of freedom on the mean field (the 
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metric of spacetime) is called backreaction, and has given rise to much speculation 
on the status of entropy and information loss in quantum gravity. A fully consis- 
tent causal formulation of the time evolution problem for quantum field theory from 
given initial conditions is absolutely necessary to discuss this wide variety of quantum 
backreaction problems. 

Apart from their backreaction on the evolution of the mean value, the fluctuating 
degrees of freedom are responsible for another important effect. In standard treat- 
ments the environment also induces a process of negative selection in the Hilbert space 
of the system: most of the quantum states become very unstable and rapidly decay 
into mixtures of relatively stable states. In other words, the environment induces 
classical behavior in the quantum system by dynamically suppressing the interfer- 
ence effects between macroscopically distinguishable states of the system. This is the 
decoherence process [8], which plays an essential role in the transition from quan- 
tum to classical behavior. As one example of this phenomenon, in a theory with 
spontaneous symmetry breaking, the coherence between the components of the wave 
function corresponding to different ground states decohere due to the coupling be- 
tween the system and the environment, so that after some characteristic decoherence 
time the state is a classical mixture (and not a quantum superposition) of macroscop- 
ically distinguishable states corresponding to different vacua. As soon as the dynam- 
ics is forced to choose between macroscopically different values (as in spontaneous 
symmetry breaking potentials) then we have the essential ingredient for macroscopic 
quantum coherence. The decoherence time scale arising from interactions with the 
fluctuating field(s) may be studied in a such a situation. The importance of deco- 
herence in the context of quantum field theory, the physics of the early universe and 
quantum cosmology has been recognized and studied in recent years [3]. Most of these 
studies focus on very simple quantum mechanical toy models and very little has been 
done in realistic quantum field theories or in closed systems generally. The techniques 
we develop in this paper are directly applicable to the study of such processes, since 
the universe as a whole is certainly a closed system. 

The first requirement for the study of non-equilibrium time evolution is a general 
initial value formulation of quantum field theory. In several earlier papers, two of 
the authors have provided the necessary Heisenberg picture formulation of the initial 
value problem in the leading order of the large N expansion [9]. It has been success- 
fully applied to the problem of pair creation in strong electric fields [10]. This leading 
order approximation in 1/N is equivalent to the Hartree-Fock mean field approxi- 
mation which has been much studied in nuclear many-body, atomic and molecular 
chemistry applications [11]. It corresonds to a Gaussian ansatz for the Schrodinger 
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wave functional. In general relativity the leading order approximation consists of 
replacing the energy- momentum source on the right side of Einstein's equations by 
its expectation value, thereby ignoring the effects of fluctuations of T M „ from its mean 
on the quasi-classical metric of spacetime. It has been clear for some time that this 
approximation is not adequate in the final stages of black hole evaporation from the 
Hawking effect, or in very early stages of the universe's expansion. However, tech- 
nical difficulties have thwarted attempts to go beyond the Gaussian ansatz for wave 
functionals. Our main contribution in this paper is a systematic technique for doing 
just that, in a way suitable for practical numerical solution on a computer. 

It is true that by suitably generalizing the Gaussian ansatz one should be able 
to evade the limitations of the mean field approximation in the Schrodinger picture 
as well. Such an approach is possible in principle but unattractive (in our view) 
principally because of the technical complications of the renormalization program 
and the loss of all connection to a covariant analysis of divergences. For example, 
in $ 4 theory, even in the lowest order approximation the Ward identities are not 
preserved by Gaussian trial wave functions without imposing the large N limit as 
well. In general, variational ansatze at the level of wave functions do not necessarily 
preserve the invariances of the underlying theory. It is the issue of renormalization on 
which most previous forays into this area founder in a web of intricate technicalities. 

In this paper we remove the restrictions of the Gaussian or mean field approx- 
imation by working in the Heisenberg picture, and making use of the Schwinger- 
Keldysh closed time path (CTP) formulation of quantum field theory [12]. It is this 
method that provides the technical means to formulate the initial value problem in 
a completely causal manner, removing the Feynman boundary conditions on Green's 
functions, which are only appropriate for (out\in) matrix elements and asymptotic 
scattering states. Otherwise, the techniques employed are completely familiar in field 
theory. Specifically, we make use of a (suitably modified) path integral representation 
for the generating function of connected Green's functions, and perform the Legendre 
transform to the effective action functional, which makes the covariance properties of 
the theory manifest. Hence, all conservation laws of the classical theory (which are not 
anomalous) are maintained explicitly. In particular, a conserved energy-momentum 
tensor for the fields in the plasma is obtained automatically by the effective action 
technique, and the renormalization program is no more difficult than in the ordinary 
covariant analysis. The variation of the effective action provides the dynamical evolu- 
tion equation(s) for the quasi-classical mean field(s) and their fluctuations, in a form 
suitable for direct numerical integration from specified initial data at t = 0. Needless 
to say these equations of motion of the (in\in) expectation values of fields are real 
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(which they would not be if the usual Feynman Green's functions were used). 

The large N expansion provides a convenient way of parameterizing the separa- 
tion into quasi-classical mean fields and their fluctuations, and at the same time, a 
systematic approximation scheme to the Schwinger-Dyson equations of the full field 
theory. The method we espouse is therefore a kind of quantum analog of truncation 
of the BBGKY hierarchy of n-particle correlation functions of classical statistical me- 
chanics at a finite order to render the system tractable to analysis [14]. Information 
loss and effective time irreversibility may enter because of this truncation [5], which 
is well motivated by the impracticality of having unlimited information about very 
high order correlations in the the initial state of any realistic system with many de- 
grees of freedom. While the large N method may not be the only approximation 
scheme that makes this separation possible, it is a gauge and renormalization group 
invariant expansion which permits non-perturbative mean fields [15]. In other words, 
the large N method goes far beyond a simple perturbative expansion in the coupling 
constant(s) which is necessarily an expansion around the vacuum field configuration, 
such as is entirely inappropriate when energy densities are high and the state is far 
from the vacuum. The large N expansion also resums and rearranges the Feynman 
perturbation series for scattering diagrams (collisions in the Boltzmann language) 
in a way that automatically includes self-energy corrections generated by the very 
same microphysical scattering and collisional processes. This points the way to a 
systematic quantal generalization of classical transport theory even in the presence 
of strong mean fields, where quantum self-energy (and all vacuum polarization or off- 
shell virtual processes) are treated on the same footing as collisional or real particle 
creation processes in the plasma. Fields which scale like TV to a positive power for 
large N can be considered strong quasi-classical mean fields in this approach. Particle 
creation effects are contained already in the leading order (i.e., N°) approximation. 
In the next to leading order (1/N) the effects of collisional, virtual and radiation 
reaction processes back on the quasi-classical fields appear for the first time. For 
almost all applications of physical interest, these effects are essential, and cannot be 
described by Gaussian wave functions. Their description requires the full power of 
the Schwinger-Keldysh initial value formulation of quantum field theory. 

In order to expose the general method in as clear as possible a manner we shall 
consider in this paper two specific field theories for definiteness, scalar A$ 4 theory 
and quantum electrodynamics (QED). Each of these is interesting in its own right. 
The A$ 4 theory may be applied to the study of the time development of phase tran- 
sitions characterized by a scalar order parameter, whether in the early universe or in 
a laboratory environment. In particular, the method is directly applicable to the chi- 
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ral phase transition and the evolution of a disoriented chiral condensate in collisions 
of relativistic heavy nuclei [17]. Electrodynamics is interesting first because QED 
is arguably the most completely verified theory of fundamental interactions known 
and an excellent testing ground for general notions of dissipation or decoherence, and 
second, for application to particle production processes in strong field astrophysical 
plasmas. Although the details of the technique will be discussed in scalar $ 4 theory 
and QED in this paper, we would like to emphasize that the technique itself is ideally 
adapted to addressing the interesting fundamental issues of quantum dissipation, de- 
coherence and time irreversibility in quantum field theories generally. The extension 
of these same methods to non-abelian gauge theories such as quantum chromodynam- 
ics (QCD), for application to the non-equilibrium time evolution of the quark-gluon 
plasma, and to gravity, for application to particle production processes in the very 
early universe and by black holes is planned in future publications. 

We are encouraged to believe that our methods will be applicable to all these 
problems by an additional technical advance in the area of renormalization. In ear- 
lier papers on the initial value formulation of quantum field theory (ours included), 
an adiabatic WKB expansion in the time variation of the frequency of oscillation 
of the fluctuating quantum modes is employed [10] [18]. The ultraviolet divergent 
contributions in the effective equations of motion appear in the first few orders of 
this asymptotic expansion and may be removed by explicit subtractions. In trying 
to extend this method beyond the lowest order mean field approximation and/or to 
spatially non-uniform mean fields we found that it quickly becomes very complicated, 
and ill-suited to practical calculations. Fortunately, we are able to dispense with this 
adiabatic expansion entirely, by the simple device of introducing an ultraviolet cutoff 
on formally divergent integrals (which one always does in practice on the computer 
in any case). Then we have only to check that rescaling the cutoff and flowing the 
coupling constant according to the standard continuum renormalization group to a 
given order in 1/N leaves physical, renormalization group invariant quantities un- 
changed. If the scale of non-equilibrium time evolution and other physics of interest 
is far from the cutoff scale this procedure is a very sensible one on physical grounds, 
and certainly it is much easier to verify that this procedure works a posteriori than it 
is to apply adiabatic type expansions beyond the lowest order in 1/N. We give some 
explicit numerical evidence of the practical feasibility of this method in this paper. 

In order to make the paper as self-contained and readable as possible, we review 
in the next two sections first the derivation of the effective action for the mean fields 
in the large N expansion for both $ 4 theory and QED, and second, the basics of the 
Schwinger-Keldysh real time CTP formalism, from the functional integral point of 
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view. In Section 4 we bring the two ingredients together and derive the causal equa- 
tions of motion for both scalar $ 4 theory and electrodynamics in the large N expan- 
sion. Section 5 is devoted to consideration of renormalization issues for the equations 
of motion, conserved currents, and energy-momentum tensors of each theory. We 
present here preliminary numerical evidence of the practicality of the renormalization 
scheme without the cumbersome adiabatic expansion employed in earlier papers. Our 
conclusions are summarized in Section 6. 



2 The Large N Expansion 

Scalar field theory with a A$ 4 self interaction is the simplest renormalizable quantum 
field theory with which to develop the techniques for separating a closed system into 
mean fields and the fluctuations about them. This model is interesting in its own 
right for the study of phase transitions characterized by a scalar order parameter, 
and for its role in the Higgs sector of the standard model of electroweak interactions. 
Let us begin by reviewing the large N expansion in the context of an TV-component 
real scalar field with the O(N) invariant Lagrangian density [15], 

L = -|$ l G- 1 [X]$ l + yx(|- M 2 ), (2.1) 

where i = 1 , . . . , N and 

G _1 [X] = -□ +X . (2.2) 
This form is equivalent to the Lagrangian density, 

v = -!(0 o $o(0°$O - - A(^) 2 (2-3) 
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with the definition of the composite field X by 



X = fi 2 + (2-4) 

since the two Lagrangians L and V differ only by a constant and a surface term. The 
quartic coupling in the Lagrangian has been taken to be X/N from the outset, rather 
than rescaling it later by 1/N as is sometimes done [16]. 

Adding independent sources to the Lagrangian and integrating over each of the 
fields X and $i defines the generating functional, 

Z[J,K] = exp(iNW[J, K]) 

= J [VX] n / [d$i] exp ji J d 4 x (£[$, X] + + NKX)} . (2.5) 



Here factors of N have been inserted which multiply both the X source K and the 
generating functional for connected Green's functions W, in anticipation of the large 
N expansion which we wish to examine. Because of the introduction of the composite 
field X, the integrations over the N scalar fields $i are Gaussian and may be performed 
explicitly: 



where the notation G[X]ij(x,y) = 6ijG[X](x,y) for the inverse of G _1 [X] has been 
used. If we take each of the N copies of the original scalar field $ to be equivalent, 
we may set each of the N components of the source Ji to be equal, and write 



in a more compact notation. In all of the following, the symbol o will denote sum- 
mation over internal indices and integration over continuous spacetime coordinates 
in the quantities on either side of it (the de Witt summation convention), whereas 
omission of the o will denote simple multiplication without summing or integrating 
over the coordinate labels. 

Once the generating functional (2.6) has been obtained by integrating over the N 
copies of the scalar field, (2.7) is obviously equivalent to simply rescaling the source 
J for a single <j> field (and the <j> field itself) by \/N. This is the sense in which the 
(j) field is quasi-classical, since its mean value is strongly enhanced with respect to 
ordinary perturbative treatments of the quantum field. A large value of the classical 
field strength implies that we are expanding about a field configuration far from the 
perturbative vacuum, which is the reason that the large N expansion is useful for 
matter under extreme conditions of high density or temperature. We should remark 
in passing as well, that scaling all N copies of the scalar field in the same way is 
not appropriate if one is interested in spontaneous breaking of the O(N) symmetry, 
in which case one component should be singled out and treated differently from the 
remaining N — 1 Goldstone fields. We shall ignore this distinction in the following, 
and focus on the O(N) symmetric case in order to simplify the presentation. 

It should be clear now why we introduced the factors of N as we did, for the 
exponent of the integrand in (2.6) contains an explicit overall multiplicative factor of 
N, and the X integration may be performed by the stationary phase method in the 
limit of large N. The stationary phase point of the integrand X S [J,K] is determined 
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(implicitly) by the relation, 
1 



K(x) + 



A 



(X(: 



-JoG{ ,x)G(x, ) o J+ -G(x,x) 



0, 



:2.s) 



and the second derivative of the exponent in (2.6) with respect to X(y) is —iND 1 
with 



x ,y) 



JoG{ ,x)G(x,y)G(y, )oJ--G(x,y)G(y,x] 



Hence the result of the Gaussian stationary phase integration over X is 
W[J,K] 



W(°) + —W {1) + — w {2) + 
N N 2 



\x a o (^ - p?) + if o X s + o G[X S ] o J + J Tr In G" 1 ^ 



(2.10) 



where X s is to be viewed as a function of the sources J and K through eq. (2.8) above, 
and order 1/N 2 terms have been dropped. By expanding the field X in a Taylor series 
about X s in the original integral in (2.6) it is straightforward to derive the corrections 
to W[J, K] appearing in (2.10) to any desired order in 1/N. Hence the approximation 
scheme is systematic and controlled, and respects all the invariances of the original 
Lagrangian L or V [16]. 

The mean field expectation values in the presence of the external sources are now 
given by the variations, 

8W N 8W 



X(c 
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" 8J{x) ' v ' 8K{x)' 

Note that the mean field X differs from the stationary phase point of the Gaussian 
integral at order 1/iV, namely 



1 6WV> 8X S ^ ( 1 

x = x ^n^x-°Jk +0 {n- 2 



(2.12) 



When no confusion with the integration variables in (2.5) is possible we shall omit 
the overline on the mean fields in the following to simplify the notation. 

The effective action functional may be defined in terms of the mean fields (2.11) 
by a Legendre transformation in the usual way, 



S[<j>,X] = W - Jo(j>- K oX, 
10 



(2.13) 



where J and K are to be regarded here as functionals of the mean fields by inverting 
eqs. (2.11). Performing the inversion to order 1/iV, we find 

J(x) = G , - 1 [X]^(a;) + -^Tr{G(x, joGo^oD}- 
1 SW™ 8X S 



o 



N 8X S " 8J{x) J=G-^<t> + ° {n 2 ) ' ^ 2 ' 14 ' ) 
When this is substituted into (2.13) and account is taken of eqs. (2.8) and (2.12), 
most of the 1/N terms cancel and we arrive at the relatively simple result, 

SIM = -^oG- 1 [X]o^+ixo(|- / x 2 ) + ^TrlnG'- 1 [X] + ^Trln J D- 1 [^X], (2.15) 
with 

ZT^XKa.iO = -\8\x,y) - 4>{x)G[X]{x,y)<l>{y) + l -G[X](x,y)G[X](y,x), (2.16) 

correct to order 1/N. By differentiating the definition of the mean fields in (2.11) 
and using the form of the generating functional (2.10) and the stationarity condition 
(2.8) again, it is easy to check that G[X] and [l/N)D[<j>,X] are precisely the lowest 
order connected two-point functions of the fluctuating quantum fields, propagating 
in background mean fields specified by (</>,X). In the following the explicit functional 
dependence of G[X] and D[<j>,X] on the mean fields will be suppressed, and we adopt 
the simpler notation G{x,y) and D{x,y) hereafter. 

From the explicit factor of 1 /N multiplying the two-point function of the X field 
and the last term in the effective action (2.15) it is clear that the fluctuations of the 
composite X field enter the discussion at one higher order of 1/N than the original 
scalar <j> field, whose fluctuations couple to the mean fields already at lowest order 
through the Tr In G~ x term in the effective action. The fact that X is an auxilliary 
field introduced into the discussion only for convenience and not a true propagating 
degree of freedom is reflected in the fact that there is no differential kinetic operator 
appearing in the expression for D~ x (unlike in the definition of G -1 ). The X field has 
no independent dynamics of its own, and its fluctuations are determined by those of 
the (j) field. The last term in (2.16), namely 

n[X] = - % -G[X]G[X] (2.17) 

is the one loop vacuum polarization of the X field due to its interaction with the 
fluctuating (j> field. It is included at the same order as the point vertex 1/A in the 
large N method because all loops of the (j> field are exhanced by a factor of N (for the 
N identical (j> fields) relative to what one would expect in an ordinary perturbative 
loop expansion. Finally we remark that 

T,{x,y) = iG{x,y)D{y,x) (2.18) 
11 



carries the interpretation of the one loop <j> self-energy, as may be seen by calculating 
the (j) inverse propagator function to order 1/iV, 

<T^X] = G" 1 + -^E + (^) . (2.19) 

The functional form of the quantum effective action (2.15) is the starting point for 
the analysis of the dynamics of the theory in the large N expansion. By differentiating 
S with respect to the mean fields we obtain their equations of motion. In the absence 
of external sources these read: 

(-□+X(xMx) + ^ Jd 4 yZ(x,y)</)(y) = 0, (2.20) 

and 

X( I ) = A^ 2 W-yG( I , I ) + ^ Jd 4 Xl J d A x 2 G{x,x 1 )Y 1 {x 1 ,x 2 )G{x,x 2 ) (2.21) 
with 

E(xi,x 2 ) = T,(x 1 ,x 2 ) - (/)(x 1 )D(x 1 ,x 2 )(/)(x2) 

= D(x 1 ,x 2 )[iG{x 1 ,x 2 )-(l>(x 1 )(l>(x 2 )] (2.22) 

The equation (2.21) for X will be recognized as just the expectation value of the 
(operator) definition of X, eq. (2.4) computed to order 1/N. 

Since the original Lagrangian (2.1) possesses the symmetry (j) — > —(j), eq. (2.20) 
is homogeneous in the <j> mean field and always admits the solution <j> = 0. In this 
form it is suitable for study of second order phase transitions. If terms breaking the 
(j) — > —(j) discrete symmetry are added (for example by retaining the linear source term 
J o 0), then first order transitions may be studied as well. Both equations for the <j> 
mean field and the X mean field contain non-local self-energy effects, and are therefore 
miegro-differential equations. In the case of the auxilliary X field the differential term 
is absent, and eq. (2.21) is an equation of constraint (or gap equation) rather than a 
true propagating equation of motion for an independent degree of freedom. 

In this derivation of the equations of motion for the mean fields from the quan- 
tum effective action (2.15) we have encountered the two-point functions G and D. 
The equations of motion should be solved concurrently with those for the two-point 
functions obtained by inverting eqs. (2.2) and (2.16) respectively. Since there is no 
unique inverse of these relations, the question of which propagator function(s) should 
be chosen presents itself. If the standard Feynman propagator functions are substi- 
tuted uncritically into eqs. (2.20) and (2.21), we find that the equations are both 
complex and acausal, in the sense that the integrations in (2.20) and (2.21) have sup- 
port in regions of spacetime that are spacelike with respect to x. The absence of a 
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well-posed causal initial value problem and the complex valuedness of the field equa- 
tions both signal that this choice of Feynman boundary conditions for the Green's 
functions is not the correct one. The reason is that we do not wish to consider off- 
diagonal (out\in) matix elements of field operators but diagonal (in\in) expectation 
values. If the field operators are Hermitian, such diagonal matrix elements must be 
real. If one's next thought is to try purely retarded propagators in the equations of 
motion (2.20) and (2.21), the result will be no better. The equations are intrinsically 
non-linear, so that no one simple choice of particular solutions to the linear equations 
G~ x o G = 1 and D~ x o D = 1 for the propagator functions will yield real, causal 
equations of motion for the mean fields. This situation is familiar in non-relativistic 
condensed matter applications, and is the reason that a more complete method for 
deriving the correct boundary conditions as well as the equations for the mean fields 
must be introduced. Having pointed out the generic non-local structure of the field 
equations following from the quantum effective action, we defer the discussion of the 
Schwinger-Keldysh CTP formalism which is precisely one such method until the next 
section, preferring first to carry out the derivation of the large N equations of motion 
in another interesting field theory, quantum electrodynamics. 

For QED with N identical charged fermion fields (flavors) we begin with the 
Lagrangian 

N r — N r 

L = ~ E 7 # x ViG-^A^i -—jtfx F ^ F>1V (2-23) 
i=i 

where anti-symmetrization with respect to the Dirac field operators ^ and ^ is un- 
derstood and 

G- 1 [A] = *y (d M ~ d„) + TA„ + im . (2.24) 
The Dirac matrices here obey 

TT + TT = 2g^ = 2 diag(-, +, +, +) (2.25) 

so that g° = — 7o is anti-hermitian and 

* = * f 7° . (2.26) 

Introducing external sources for the gauge potential and Dirac fields, we define the 
generating functional 

Z[J,r,,rj] = exp(iNW[J,K]) = J [VA^] 1 J[ J[dVi][d%] exp ji J d 4 x L[A, *]| 

x exp \iNJ o A -rjo $ - * o 77} , (2.27) 

where the prime on the gauge field integration measure denotes that we should inte- 
grate only over distinct gauge invariant configurations (or equivalently, fix the gauge). 
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Performing the Gaussian integration over the anti-commuting Dirac fields, and rescal- 
ing the Grassman valued sources 77 — > \fNrj so that we can drop the sums over 
i = 1, . . . , N as in the previous scalar case, we obtain 

Z[J,V,V] = JiVA^'exp^Njd'x^A^g^D-d^A^ 

x exp [NTv In G" 1 - Nrj o G[A] o 77 + iNJ o A) . (2.28) 

As in the previous discussion we have defined the sources and coupling with the 
correct powers of N to justify performing the remaining functional integration over 
the electromagnetic potential by the stationary phase method. The stationary phase 
value A 8 ^ [J, 77,77] is fixed by 

^{g^U- d^d") A s v {x) = itr {G[A s ](x, x)T}+irjoG[A s ]( , x)TG[A']{x, )o V -J»(x) 

(2.29) 

where tr denotes the Dirac matrix trace only (without integration over spacetime 
coordinates). The second derivative of the exponent in (2.28) at its stationary point 
is —iNTrD -1 where 

D-'iJ^M^vr = -^(g^a-d^)s\x,y) 

- ztr {G(y, x)T G(x, y)T + 2rjoG(x)TG(x,y)TG(y, ) o 77} (2.30) 

and A = A S [J, 77,77] from the eq. (2.29) and symmetrization with respect to inter- 
change of the pair of spacetime labels (x,/x) with {y,v) is understood. Thus, the 
result of the stationary phase evaluation of (2.27) and (2.28) in the QED case is 

W[J, K] = -A s o d' 1 o A s + J o A s + irf o G[A S ] o 77 

-zTrlnG- 1 ^] + ^Tr In IT 1 [A s ], (2.31) 

where order 1/N 2 terms have been dropped, and 

d-\x,yr = -^(9^0 - d^d v )6\x,y) (2.32) 

is the differential operator from the classical action. 

The mean fields and quantum effective action are defined now by the analogs of 
eqs. (2.11), and (2.13), after solving for the sources 7,77,77 in terms of the mean fields. 
Omitting the details which are quite analogous to the scalar case, the result of this 
Legendre transformation is simply 

S[A] = -^Aod' 1 o A-iTrlnG-^A] + ^Tr In D' 1 [A], (2.33) 
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in the case of zero mean value for the Dirac field. The photon inverse propagator in 
the last term is given by 



D-'lA^x^r = (d- 1 + U[A])(x, y r (2.34) 

with 

Il[A](x,yy» = -itT{TG[A](x,y)TG[A](y,x)} (2.35) 

the polarization tensor in the presence of the mean potential A. The inverse prop- 
agator cannot be inverted without fixing a gauge, which may be done by a variety 
of standard methods. In the case of a non-abelian gauge symmetry, the gauge fixing 
introduces ghosts which will also contribute to the quantum effective action at order 
I/TV, and whose contribution is essential to obtaining gauge invariant results. In the 
abelian QED case the ghosts are independent of the mean potential (i.e., they de- 
couple) and therefore may be omitted from the effective action (2.33). Of course, the 
non-abelian case is of great interest for describing the non-equilibrium evolution of 
the quark gluon plasma from first principles of QCD, and will be discussed in detail 
in a future publication. 

The integro-differential equation for the mean potential is obtained by varying the 
effective action. In this variation will appear 



+itr [l a G(x 1 ,x 2 )l^G(x 2 ,x)TG(x,x 1 )] , (2.36) 



so that the equations of motion for the mean potential read 



2 



d v F^(x) = —ie tr {7 M G[A](x, x)} H — Tr < D[A] o 



m[A] 



2N I 8A„ 



x 



-% 



e 2 tr{7 M G(x,x)} 



+ J d 4 Xl J d 4 x 2 tr {l»G(x, x 1 )E(x 1 , x 2 )G(x 2 , x)} (2.37) 



wi 



ith 



E(xi, x 2 ) = iTG(x 1 ,x 2 )rD vll (x 2 , xi) (2.38) 

the fermion self-energy. The current expectation value to order 1/N is represented 
pictorially by the graphs in Fig. 1. 

To leading order in 1/N the mean field equations are just the semi-classical 
Maxwell equations, obtained by replacing the electric current operator of the Dirac 
field by its expectation value. This leading order semi-classical equation already con- 
tains the dynamical reaction of e + e~ pairs created by a non-zero electric field (the 
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Figure 1: First two graphs in the 1/N expansion contributing to the induced current 
that governs the back reaction on the electric field. 

Schwinger mechanism) back on the electric field itself, and has been studied in previ- 
ous publications [10]. However, at leading order in 1/N the created pairs can interact 
only through the mean field, not directly with each other. The order 1/N term with 
the fermion self-energy E contains the quantum Compton scattering, bremmstrahlung 
and Coulomb interaction effects of these particles on each other and their backreaction 
on the self-consistent mean field. Clearly these processes are essential mechanisms for 
the approach to equilibrium and must be included in any realistic transport theory of 
a relativistic e + e~ plasma. Higher order processes (such as multiple scattering) can 
be included by working to higher orders in the 1/N expansion in a straightforward 
way. 

Having derived the integro-differential equations for the mean fields in both scalar 
A$ 4 theory and electrodynamics to order 1/N, we turn now to the Schwinger-Keldysh 
closed time path formulation of the effective action, in order to determine the correct 
propagator functions needed to obtain a causal (and real) solution to these equations. 

3 The Schwinger-Keldysh Closed Time Path For- 
malism 

The conventional path integral formalism used freely in the preceeding section defines 
transition elements between states at one time, t (usually taken to be in the infinite 
past) to states at another time t' (in the distant future). If the class of paths is 
restricted to be the vacuum configuration at both of its endpoints, then the two 
states are the \in) and (out\ vacuum states of scattering theory respectively. The 
functional Z[J,K] of eq. (2.5) is the transition matrix element 

Z[J, K](t, t') = {out, t'\in, t) JiK (3.1) 

in the presence of the external sources J and K . 

By varying with respect to the external sources we obtain matrix elements of the 
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Heisenberg field operators between the \in) and (out\ states. For this reason we may 
refer to the conventional formulation of the generating functional Z as the "in-out" 
formalism. The time-ordered Green's functions obtained in this way necessarily obey 
Feynman boundary conditions, and these are the appropriate ones for the calculation 
of transition probabilities and cross sections between the \in) and (out\ states. On 
the other hand the off-diagonal transition matrix elements of the in-out formalism 
are completely inappropriate if what we wish to consider is the time evolution of 
physical observables from a given set of initial conditions. As we have remarked the 
in-out matrix elements are neither real, nor are their equations of motion causal at 
first order in l/N, where direct self interactions between the fields appear for the first 
time. What we require is a generating functional for diagonal matrix elements of field 
operators with a corresponding modification of the Feynman boundary conditions 
on Green's functions to ensure causal time evolution. This "in-in" formalism was 
developed more than thirty years ago by Schwinger and later by Keldysh, and is 
called the closed time path (CTP) method [12]. 

The basic idea of the CTP formalism is to take a diagonal matrix element of the 
system at a given time t = and insert a complete set of states into this matrix 
element at a different (later) time t' . In this way one can express the original fixed 
time matrix element as a product of transition matrix elements from to t' and the 
time reversed (complex conjugate) matrix element from t' to 0. Since each term in this 
product is a transition matrix element of the usual or time reversed kind, standard 
path integral representations for each may be introduced. If the same external source 
operates in the forward evolution as the backward one, then the two matrix elements 
are precisely complex conjugates of each other, all dependence on the source drops 
out and nothing has been gained. However, if the forward time evolution takes 
place in the presence of one source J + but the reversed time evolution takes place 
in the presence of a different source J_ , then the resulting functional is precisely the 
generating functional we seek. Indeed (setting K = and N = 1 here for simplicity), 



Z m [J + ,J_] = J[V^](in\^)j_ (^\in) J+ 

= J[V^](in\T*exp -i J dtd 3 xJ_ 



x(q,t'\T exp 



/' 

Jo 



dtd 3 xJ + (x)(/)(x 



x)(j>(x) 
in) 



(3.2) 



so that, for example, 

6W m [J+,J- 



8J+(x) 



6W m [J + ,J. 



J+=J_0 



6J-(x) 
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j + =j_ 



(in\(j)(x)\in) 



(3.3) 



© 



Figure 2: Complex time contour C for the closed time path propagators. 



is a true field expectation value in the given time-independent Heisenberg state \in). 
Since the time ordering in eq. (3.2) is forward (denoted by T) along the time path 
from to t' in the second transition matrix element, but backward (denoted by T*) 
along the path from t' to in the first matrix element, this generating functional 
receives the name of the closed time path generating functional. If we deform the 
backward and forward directed segments of the path slightly in opposite directions 
in the complex t plane, the symbol 7c may be introduced for path ordering along the 
full closed time contour, C depicted in Fig. 2. 

This deformation of the path corresponds precisely to opposite ie prescriptions 
along the forward and backward directed segments, which we shall denote by C± 
respectively in the following. 

The doubling of sources, fields and integration contours in the CTP formalism 
may seem artificial, but in fact it appears naturally as soon as one discusses the time 
evolution not of states in Hilbert space but of density matrices. Then it is clear that 
whereas | ) ket states evolve with Hamiltonian H , the conjugate ( | bra states evolve 
with —H, and the evolution of the density matrix requires both. Hence a doubling 
of all sources and fields in the functional integral representation of its time evolution 
kernel is necessary. Indeed, it is easy to generalize the functional in (3.2) to the case 
of an arbitrary initial density matrix p, by defining 



Z[J+,J-,p] = Trip IT* exp 



ft' 

. / dtd xJ_(x)</)(x) 
Jo 



X T exp 



l-t' 

i / dtd xJ + (x)</)(x) 
Jo 



j ' [V<p][V<p'][VTP] (<p\p\<p' ) (<p'\T*exp -i j* dtd 3 xJ_(x)</)( 



X (ip\Texp 



ft' 

i / dtd xJ + (x)</)(x) 
Jo 



\<P) ■ 



(3.4) 
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Variations of this generating function will yield Green's functions in the state specified 
by the initial density matrix, i.e., expressions of the form, 

Tr{^(x 1 )0(x 2 )0(x 3 )...}. (3.5) 

Introducing the path integral representation for each transition matrix element in 
eq. (3.4) results in the expression, 



exp 



poo 

i \ dtd 3 x( L[(j) + ] - L[(j)_] + J + (j> + - J_(j>. 
Jo 



(3.6) 



where L is the classical Lagrangian functional, and we have taken the arbitrary future 
time at which the time path closes t' — > oo. 

The double path integral over the fields and in (3.6) suggests that we 
introduce a two component contravariant vector of field variables by 

r= ^ + j ; a = 1,2 (3.7) 

with a corresponding two component source vector, 

J a = ^ + j ; a =1,2 . (3.8) 

Because of the minus signs in the exponent of (3.6), it is necessary to raise and lower 
indices in this vector space with a 2 X 2 matrix with indefinite signature, namely 

c ab = diag (+1,-1) = c ab (3.9) 

so that, for example 

J a c ab $ b = J-<f>- . (3.10) 

These definitions imply that the correlation functions of the theory will exhibit a 
matrix structure in the 2x2 space. For instance, the matrix of connected two point 
functions in the CTP space is 

S 2 W 



G ab (x,y) 



(3.11) 

j=o 



6J a (x)6J b (y) 

Explicitly, the components of this 2x2 matrix are 
G 21 {x,y) = G > (x,y)= zTrV $(x)$(2/)} con , 

G 12 {x,y) = G < (x,y) = ±iTr{p$(y)$(x)} cm , (3.12) 
G n (x,2/) = iTr{pT[*(x)*(y)]\ = 6(x,y)G > (x,y) + 6{y , x)G < (x,y), 

< ■> con 

G 22 (x,y) = iTr{ /> T*[$(x)$(y)]} = 6(y, x)G > (x, y) + 6{x, y)G < (x, y), 

< con 
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where the ± refers to Bose or Fermi fields respectively, and $ = $t for bosons but 
\|r = ^rt-yO £ Qr j)j rac fermions. Notice that G > (x,y) = G < (x,y) if x and y are spacelike, 
and therefore 



G (x,y) = G (x,y) for (x, y) spacelike. (3.13) 

and in particular on the spacelike surface t x = t y . This will be fully consistent with 
the expressions for time ordered and anti-time ordered propagator functions in (3.13), 
provided we extend the usual definition of the discontinuous 6 function so that 



so that in particular, 



6(t,t') + 6(t',t) = l for all M' 



(3.14) 



(3.15) 



The 2x2 matrix notation has been discussed extensively in the literature [12]. 
However, the development of the CTP formalism is cleaner, both conceptually and 
notationally, by returning to the definition of the generating functional (3.4), and 
using the composition rule for transition amplitudes along the closed time contour in 
the complex plane. Then we may dispense with the 2x2 matrix notation altogether, 
and write simply 



J [Dip] (<p'\T* exp —i J 



dtd 3 x J_ 



(ip'\T c exp 



r poo 

)cj){x) \ti)){ij)\T exp i / dtd 3 xJ + {x)cj){ 
J L Jo 

J dtd 3 x J(x)(j)(x) \ip) 



W) 



(3.16) 



so that (3.4) may be rewritten more concisely in the CTP complex path ordered form, 

Zc[J,p] = Tr j^^Tcexp i J dtd 3 xJ(x)(j)(x) 

= J [Dip 1 ] J [Dip 2 ] (c/H^ 2 ) J* [D<f>] exp % J dtd 3 x (L[(j>] + J</>) 

(3.17) 

The advantage of this form is that it is identical in structure to the usual expression 
for the generating functional in the more familiar in-out formalism, with the only 
difference of path ordering according to the complex time contour C replacing the 
ordinary time ordering prescription along only C + . Hence, all the functional formalism 
of the previous section may be taken over line for line, with only this modification of 
complex path ordering in the time integrations. For example, the propagator function 
becomes 



G(x,y) = 9 c (t x ,t y )G > (x,y) + (9 c (t y , t x )G'<(x, y) 
= 6 c (t x ,t y )G 21 (x,y) + 6 c (t y ,t x )G 12 (x,y] 
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(3.18) 



where 8c is the CTP complex contour ordered theta function defined by 



e c (t,f) 



6(t,t') for t,t' both on C+ 

6{t'A) for t,t' both on C_ 

y ' J ' (3.19) 

1 for t on C- , t' on C + 

for t on C + , t' on C_ 



With this definition of G{x,y) on the closed time contour, the Feynman rules are the 
ordinary ones, and matrix indices are not required. In integrating over the second 
half of the contour C_ we have only to remember to multiply by an overall negative 
sign to take account of the opposite direction of integration, according to the rule, 

/* / i oo / i oo 

dt= dt- dt . (3.20) 
Jc Jo c + Jo C- 

A second simplification is possible in the form of the generating functional of 
(3.17), if we recognize that it is always possible to express the matrix elements of the 
density matrix as an exponential of a polynomial in the fields [5]: 

(c/|/V ) = exp R + J d 3 xR a (x)ip a (x) + J d 3 xd 3 yR ab (x, y)ip a (x)ip\y) + ... . 

(3.21) 

Since any density matrix can be expressed in this form, there is no loss of generality 
involved in expressing p as an exponential. If we add this exponent to that of the 
action in (3.17), and integrate over the two endpoints of the closed time path if 1 
and if 2 , then the only effect of the non-trivial density matrix p is to introduce source 
terms into the path integral for Zc[J,p\ with support only at the endpoints. This 
means that the density matrix can only influence the boundary conditions on the 
path integral at t = 0, where the various coefficient functions R a , R a b-, etc. have the 
simple interpretations of initial conditions on the one-point (mean field), two-point 
(propagator), functions etc. It is clear that the equations of motion for t ^ are 
not influenced by the presence of these terms at t = 0. In the special case that 
the initial density matrix describes a thermal state, p@ = exp{— f3H} then the trace 
over pp may be represented as an additional functional integration over fields along 
the purely imaginary contour from t = —if3 to t = traversed before C_ in Fig. 2. 
In this way the Feynman rules for real time thermal Green's functions are obtained 
[13]. Since we consider general nonequilibrium initial conditions here we have only the 
general expression for the initial p above and no contour along the negative imaginary 
axis in Fig. 2. 

To summarize, we may take over all the results of the previous section on the 
generating functionals, effective actions, and equations of motion of scalar $ 4 theory 
and QED, provided only that we 
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1. substitute the CTP path ordered Green's function(s) (3.18) for the ordinary 
Feynman propagators in internal lines; 

2. integrate over the full closed time contour, C, according to (3.20); and 

3. satisfy the conditions at t = corresponding to the initial density matrix p. 

4 Causal Evolution Equations in <I> 4 Theory and 
QED 

To show that these three modifications of ordinary in-out formalism do lead to a well- 
posed initial value problem for quantum field theory with real and causal equations 
of motion for the mean fields and their Green's functions, let us reconsider eqs. (2.20) 
and (2.21), Section 2. First, the self-energy has the CTP ordered form, 

T,(x,y) = i6 c (t x ,t y )G > (x,y)D > (x,y) + i6 c (t y ,t x )G < (x,y)D < (x,y) (4.1) 

since both G{x,y) and D{x,y) have this form separately. Using the rule (3.20) for 
the time contour integration, we have 

dt y Ti(x,y)^(y) = i dt y G > {x,y)D > {x,y)cj){y) + i dt y G < {x,y)D < {x,y)cj){y) 
Jc Jo Jt x 

poo 

-i / dt y G < (x,y)D < (x,y)(j)(y) 
Jo 

= i dt y [G > (x,y)D > (x,y) - G < (x,y)D < (x,y)] (j)(y) 
Jo 

= -2 t dty lm[G > (x,y)D > (x,y)]<l>(y) (4.2) 
Jo 

The minus sign in the second line is because t y is on the second branch of the contour, 
and the fact that 

[G>(x,2/)r = -g<(*,io 

[D>(x,y)]* = -D<(x,y) (4.3) 

for real boson fields has been used. Therefore, the equation of motion for the mean 
field, (2.20), is 

(-□+X(xMx)--| fj dt y d 3 y lm[G > (x,y)D > (x,y)]</>(y) = . (4.4) 

This equation is now explicitly both real and causal. Furthermore, it is not difficult 
to see that the cancellation of the acausal parts of the integration (when t y > t x here) 
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is completely general in the CTP method, since every internal vertex requires a time 
integration over the full path C, and equal and opposite contributions will always 
come from the portions of the integration on C + and C_ with t y > t x . One simply 
decomposes the time integration for each internal vertex t{ into three segments, viz., 

(i) < U < t, on C + 

(ii) t < ti < oo, on C + 

(m) < U < oo, on C_ (4.5) 

uses the definitions (3.18), (3.19), (3.20), and collects the non-cancelling terms. In 
this way we find that the 1/N term appearing in the eq. (2.21) for the X mean field 
reduces to 

% J dt! J dt 2 G{t, t2)G{t 1 ,t)t(t 1 ,t 2 ) = i J* dt t J* dt 2 {G > {t, t 2 )G < {t 1 ,t)f,{t 1 ,t 2 ) 

- G<(t,t 2 )G' < (t 1 ,t)S < (t 1 ,t 2 ) 

- G>(M2)G>(*i,0S>(*i,*2) 

- G < (t,t 2 )G > (t 1 ,t)f:*(t 1 ,t 2 )}, (4.6) 

where the last term is anti-time ordered since t\ and t 2 are both on C_ , and all spatial 
dependences have been suppressed. By now expanding in the possible orderings of t\ 
and t 2 , interchanging integration variables, and using 

[G>,<{x,y)\* = -G >t< (y,x) 

[D >t< (x,y)]* = -D >t< {y,x) (4.7) 
which is true for real or complex Bose fields, we can write (4.6) in the form, 

2 Jm y V p <ft 2 [G>(*i,*) - G<(*i,*)] [G>(t,t 2 )S > (t 1 ,t 2 ) - G<(t,t 2 )S < (t 1 ,t 2 )] 

(4.8) 

which is manifestly real and causal. It is easy to see that this result is identical to 
that obtained in the matrix notation from 



; f dt x { 

Jo Jo 



1 ' dt 2 G 1 \t,t 2 )G a \t 1 ,t)Y, ab (t 1 ,t 2 ) (4.9) 



provided indices are raised and lowered with the indefinite metric c a b defined in eq. 
(3.9). Since the causality and reality of the equations are assured on general grounds, 
the result of any calculation is obtained most rapidly by proceeding to this matrix 
form directly. 

The effective Feynman rules for the CTP method in the large N expansion may be 
summarized then as follows. Derive the effective action and equations of motion for 

23 



the mean fields and propagator functions to a given order in 1 /N in the usual in-out 
formalism, as in the previous section. In the corresponding graphical representation 
of these equations, replace the propagators with the 2x2 matrix propagator of the 
CTP method, raising and lowering all contracted indices at internal vertices with 
the metric c a {,, and fixing the "external" vertex of the mean field(s) to be of type 1. 
After all time orderings have been taken into account only the Wightman functions 
such as G> = G 21 or G < = G 12 will appear in the equations, which can be solved 
for concurrently with the corresponding mean fields in a well-posed real and causal 
intitial value problem. 

For the $ 4 theory we have gone through the CTP procedure in some detail using 
the above manipulations. The final result for the equation of the <j> mean field has 
been given already by (4.4). The corresponding causal X equation (2.21) is 

X(x) = fi 2 + ^<j> 2 (x)-^[G > (x,x) + G < (x,x)] 

\ rt yti 

+ — Xm / dtid 3 Xi / dt 2 d 3 x 2 [G > (x, X\) — G < (x, x\)] 
N Jo Jo 

x [S<(xi, x 2 )G > (x 2 , x) - S>(xi, x 2 ) - G<(x 2 , x)] , (4.10) 

with E given by eq. (2.22) of the previous section. 

Concurrently with the mean field equations we must solve for the two-point func- 
tions G and D. If both have the causal form of (3.18) in the path ordered notation, 
then the Wightman functions G >t< satisfy the homogeneous equations, 

(-□ +X(x))G > ,<(x,x') = , (4.11) 

together with the initial conditions following from the initial density matrix p as 
specified by eqs. (3.13). Because of the canonical equal time commutation relations 
the initial conditions on the first derivatives of G >t< must satisfy the constraint, 

^ [G^y) - G^y)]^ = 8\x-y) (4.12) 

which guarantees that the inhomogeneous equation G~ x o G = 1 is satisfied by (3.18). 
There is no corresponding derivative condition for the D propagator function since the 
operator D~ x of (2.16) does not contain any time derivatives. Instead, the function 
D satisfies the inhomogeneous integral equation, 

8\x,y) = -jD{x,y) - d l x 1 fl(x,x 1 )D(x 1 ,y) (4.13) 

with 

U(x,y) = IL(x,y) + <f>(x)G(x,y)<f>(y) 

= G(x,y)[- l -G(y,x) + ( l ) (x) ( l ) (y)} ■ (4-14) 
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The only way to satisfy the inhomogeneous equation (4.13) is for the D y < functions 
themselves to contain a 8 function term. This is simply a consequence of the fact 
that the X field defined by eq. (2.4) is an auxilliary field, purely constrained in terms 
of (j) by its definition with no conjugate momentum or independent dynamics of its 
own. Hence we seek a solution of (4.13) of the form, 

D >t< {x,y) = -^\x,y) + D >t< {x,y) (4.15) 
with D y < smooth functions satisfying 

D<,>(x,y) = X 2 Il <j> (x,y) - A / d 4 x 1 [Il > (x,x 1 ) -Il > (x,x 1 )]D <t> (x 1 ,x) 

Jo 

+ A f ty d 4 x 1 Il <j> (x,x 1 )[D > (x 1 ,y) -D < (x 1 ,y)] . (4.16) 
Jo 

Evidently the functions D y < are determined completely by the initial conditions and 
causal evolution of the (j> field and its propagator functions G> i<) consistent with X 
being a fully constrained field with no independent dynamics. 

This completes the derivation of the causal equations of motion for the initial 
value problem of nonequilibrium $ 4 theory to order 1/N. The derivation of the 
corresponding equations for nonequilibrium electrodynamics proceeds in an exactly 
analogous manner with the only difference that the relation 

[G >t< {*,yW = -7 G > , < ( 2/) x)7 (4.17) 

for complex Dirac fields replaces the first member of (4.7). The result is that the 
causal Maxwell equations of motion take the form, 

d v F^{x) = (j»(x)) = --^-tr{T[G>(x,x) + G < (x,x)}} 

2e 2 t l /"'i f 

+ / dt x d z x x \ dt 2 d 3 X2lmtr\l^ [G>(x,Xi) - G<(x,Xi)] 

N Jo Jo I 

x [E < (x 1 ,x 2 )G > (x 2 ,x) - E > (x 1 ,x 2 )G < (x 2 ,x)] | . (4-18) 

The Wightman functions for the Dirac field satisfy 

{Td» - iTAp + m) = (7"V M + m)G >i< {x,y) = 0, (4.19) 

together with the initial conditions implied by the first two members of eqs. (3.13), 
which satisfy the canonical equal time ant i commutator condition, 

- i [G>(x >y ) - G<(x,2/)] ts=ty = Tr {,[¥(*), *(2/)] + } ts=ty = " V) > ( 4 - 20 ) 

appropriate for Fermi-Dirac statistics. 
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For a complete initial value problem to order 1/iV, one needs also the two-point 
function of the Maxwell field obtained by inverting (2.34) subject to some gauge 
condition. The simplest way to impose the gauge condition is to treat it in a way 
similar to the constraint of the A$ 4 theory, i.e., we write 

D nv(x,y) = e 2 d llv {x,y) + D^,{x,y) (4-21) 

with d^v the inverse of the differential operator (2.32) of the free action in a definite 
gauge. The gauge fixing can be performed at the level of the free photon propagator 
once and for all, independently of the dynamical time evolution problem, and the 
non-trivial time evolution is contained entirely in D^, which obeys 

( g^U - d»d X ) D Xv (x,y) = e 2 J d 4 Xl W\x, Xl ) (e 2 d Xv { Xl ,y) + D Xv ( Xl ,y)) . 

C (4.22) 
Writing the propagator in terms of its time ordered structure along the contour C as 
in eq. (3.18) we obtain 

(g^U-d»d X )Dll < { X ,y) = 

-e 2 fj d 4 Xl H^ i< (x,x 1 ) le 2 d>„( Xl ,y) + D> v { Xl ,y) - fd^x^y) - 5<„(xi,y)] 
+e 2 jf'Vx! [nf^xO-n^x,^)] {e 2 dt < ( Xl ,y) + D>f( Xl ,y)} , (4.23) 

which is explicitly causal. We remark that a particularly useful gauge choice for 
practical implementation of the initial value problem on a computer is the Coulomb 
gauge, which has the advantage of clearly isolating the physical transverse modes of 
the photon and allowing the longitudinal and gauge modes to be eliminated from the 
evolution problem, thereby making most efficient use of computer memory. In the 
Coulomb gauge the propagator and each of its pieces satisfies 

d m D mv = d m d mv = d m D mv = . (4.24) 

The Coulomb gauge condition substituted into (4.23) allows one to solve for the time 
components of the propagator D tv explicitly in terms of the other components since 
there are no propagating timelike photons in this physical radiation gauge. 

In cases of special symmetry, such as spatially homogeneous mean fields, con- 
siderable simplification of the equations of motion occur. In the scalar field case let 
X = X(t) and <j> = <j>(t) be functions only of time. Then the Wightman functions may 
be expressed as the Fourier transform, 

G>(t, x; t', x') = -G%(t, x; t', x') = J [dk]^*-*'^^, t'; k) (4.25) 
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with 

W = p^I (4-26) 

and G> a function only of k = |k|. Introducing the Fourier mode functions fk(t) 
satisfying 

^ " /*(*) = (4-27) 



d2 ,2 



and normalized by the Wronskian condition, 

we may express the Wightman function in the form, 

G>(t, t'; k) = zf k (t)f* k (t')(N(k) + 1) + ifl{t)f k {t')N{k) + 2z^e (/ fc «/ fc (0^)) 

(4.29) 

where N(k) and ^(A;) carry the interpretation of particle number and correlated pair 
density in the general spatially homogeneous initial state: 



Tr(^4«k') = (2vr) 3 £ 3 (k - k')N(k) 

TrUi k a k ,) = (2vr) 3 £ 3 (k + k')F(k) . (4.30) 



By likewise Fourier transforming the D propagator, the field equations for the spa- 
tially homogeneous case may be expressed in the final form, 



^ + X{t)\ 4{t) = ^ J*dt' J o °° kHklmiG^t'-^D^t'-k)]^') , (4.31) 



A A f°° 

X(t) = + -<j)\t) + — / k 2 dk 1m G>(t, t; k) + (4.32) 

I 47T^ JO 

!l dtl So dt2 So" k2dkne [G>i t , t i\ k )] Im p>(*i,*2; Aj)G>(* 2 ,*; A)] , 



with 



_ _ r t' 

= A 2 n>(t,t'; Jfe) + 2A / dt^^^i, k)1Ze Dy^^'; k) 

Jo 

-2A^ [fte n>(Mi;&)] Dy^t^t'-k) (4.33) 

n>(M'^) = ^)G > (t > ^A ! )^(0-^/[dq]G > (t > ^ g )G > (t' > t;|k-q|) 
£>(M';&) = -^(0^>(*,* , ;A : )^ , ) + iy[dq]G > (* > * , ; 9 )5' > (* > * , ;|k-q|) 

(4.34) 
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and eqs. (4.27)-(4.29) above. All the equations have been expressed in terms of the 
G> and Dy functions by using the relations (4.3) and (4.7). 

The result of our derivation is a well-defined initial value problem posed by the 
closed set of causal equations of motion and initial data, (4.27)-(4.34), corresponding 
to the mean fields and their fluctuations evolving forward in time from a specified 
spatially homogeneous initial density matrix. These equations are now suitable for 
numerical solution, provided the ultraviolet divergences in the momentum integrations 
are absorbed into renormalized parameters of the theory in the usual way. It is to 
this technical issue of the renormalization procedure that we turn next. 

5 Renormalization and Energy- Momentum Ten- 
sor 

To lowest order in 1/N (which includes the time dependent Hartree-Fock or Gaussian 
approximation) a convenient approach to the identification and removal of ultraviolet 
divergences in the evolution equations is the adiabatic method. This method is based 
on the fact that the large momentum behavior of the mode functions satisfying the 
differential eq. (4.27) is determined by an asymptotic expansion in derivatives of the 
time dependent frequency ^k 2 + X(t). The divergences in the currents appearing in 
the mean field equations for X(t) or A^(t) are contained in the first few terms of 
this adiabatic expansion for the mode functions, and coincide with the divergences 
of the manifestly covariant in-out formalism. Hence they are removed by the same 
countertems as in the covariant approach. 

At next order in 1/N the structure of the integro-differential equations for the 
mean fields is considerably more complicated, and the adiabatic method for identify- 
ing and explicitly removing the ultraviolet divergences appears to be quite unwieldy. 
In solving the evolution equations numerically, an ultraviolet (and infrared) cut-off 
is always present. There is no real need to remove the explicit cut-off dependence 
appearing in intermediate quantities, provided only that physical results are cut-off 
independent in the end. If all bare quantities are taken to depend on the explicit 
ultraviolet momentum cut-off and the renormalized parameters of the theory in the 
same way as in the usual covariant treatment, then the results for the physical time 
evolution of the system should be insensitive to the the cut-off in the final analysis. 
This insensitivity of the solution to the cut-off may be checked empirically by chang- 
ing the cut-off and evolving the equations again with the same initial data. That is, 
we increase the cut-off and rescale the bare parameters keeping the renormalized pa- 
rameters fixed until the evolution is insensitive to the cut-off. On the order of several 
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thousand field modes are typically required to approach this regime of insensitivity 
to the cut-off. 

To illustrate the practicality of this method consider the simplified case of spatially 
homogeneous mean fields in the lowest order of the 1/N expansion. The scalar $ 4 
equations read in this case simply, 



d 2 

dT 2+X ^ 



m = o 

A 



X(i) = / . 2 + ^) + -4lm/ k 2 dk 0^1,1; k), (5.1) 
2 vk 1 Jo 

with G>(t,t; k) given in terms of the mode functions fk(t) by eqs. (4.27)-(4.29), and 
the dependence of the bare coupling upon the momentum cut-off is recorded explicitly. 
The X equation contains a quadratic divergence in the momentum integration as 
the cut-off A — > oo, which is independent of time and must be compensated by 
a counterterm in the bare mass parameter \i 2 . This quadratic divergence may be 
removed by the simple device of expressing X(t) in terms of its finite initial condition 
at t = 0, i.e., 

X(t) = X(0) + y (<f>\t) - 2 (O)) + ^ f k 2 dk lm {G>(t, t; k) - G>(0, 0; *)} • 

(5.2) 

The integration over k still leads to a logarithmic dependence on the cut-off for large A, 
but this is precisely compensated by the A dependence of the bare coupling according 
to the usual renormalization framework, viz., 

A A = Z^(A,m)X R (m 2 ) (5.3) 

with 



Z x (A,m) = 1 - — -X R (m 2 )ln (-) 
loir 2 \mj 



l + -^A A ln(-)l 1 (5.4) 
loir 2 \mJl 

and A#(m 2 ) the renormalized quartic coupling defined at some finite mass scale m 2 . 
Indeed, by dividing both sides of (5.2) by A A and using eqs. (5.3) and (5.4), it is 
straightforward to verify that the logarithmic dependence on A of the integral in 
(5.2) is cancelled by the logarithm in (5.4), so that the resulting equation for X(t) is 
independent of A for A large (provided Z\ > 0). Since there is no <j> (wavefunction) 
renormalization at lowest order in 1 /N these simple steps are all that are required to 
arrive at cut-off independent evolution equations for the scalar theory at this order. 
Notice, in particular, that no adiabatic expansion of mode functions to isolate diver- 
gences is necessary in this approach. Such an expansion is useful only for verifying 
explicitly the cancellation of the cut-off dependence in (5.2) which must occur in any 
case. 
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In QED the Maxwell equation for the homogeneous mean electric field in the 
gauge A = and A{ = 6i z A(t) at lowest order in 1/N becomes simply 

J2 ' 2 

—A{t) = e\Tr{pf{t)) = --^ y[dk]tr{[G > (t > t;k) + G<(t > t;k)]7*} . (5.5) 

The Wightman functions G >t< may be expressed in terms of spinor mode functions 
obeying 



through 



G>,<(M';k) 



i— a ■ (k — A) — Bra 

dt y ' H 



iE.^k.W^k.^') 








(5.6) 



+ 



( -*£.«k.(0 T, k.(* / ) 

S 

+z £ (F(ks)u ks (t)v ks (t') + F*(ks)v ks (t)u ks (t')) (5.7) 

s 

respectively, where N and F are the mean number of particle pair and correlation 
densities in the initial state, and 5 = 1,2 labels the spin. By assumption, N and F 
fall off faster than |k|~ 4 for large |k|, so that the only divergence on the right side of 
eq. (5.5) comes from the N = F = (vacuum) contribution, viz., 

eA<ob'*(Olo) = y /[dk]E {^.(0^.(0-^.(0^.(0} • (5-8) 

A naive cubic divergence in this expression is cancelled between the two charge conju- 
gated spinors u and v. There remains only the logarithmic cut-off dependence related 
to charge renormalization. Let 



(5.9) 



Wl 



ith 



Z e (A,m) = 1 



6-7T 



1 eifmW^ 



2 c iA 



raj 



1 + 



1 



6-7T 



elln 



2 C A 



ray 



(5.10) 



which is the usual covariant charge renormalization of QED to this order. By dividing 
both sides of the Maxwell equation of motion, (5.5) by e\ and using the relations (5.9) 
and (5.10) it is straightforward to check that the logarithmic cut-off dependence of the 
current expectation value is precisely cancelled by the logarithmic A dependence of 
(5.10). Numerical results demonstrating the cut-off independence of the final result 
for the time evolution are presented in Figs. 3. The adiabatic expansion of mode 
functions previously employed in Refs. [9] [10] [18] is therefore quite unnecessary. 
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Figure 3: (a) Unrenormalized time evolution of the electric field in lowest order in 
1/N for fixed initial conditions and fixed bare charge (e 2 = 50) but different values 
of the cutoff. The electric field is scaled by the critical field (E — > eE/m 2 ), and the 
time and cutoff are scaled by the mass (the cutoff is given in units of m while time is 
measured in units of m -1 ). The solid, dashed, and dotted lines correspond to cutoffs 
of 64, 48, and 32, respectively. 
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Figure 3: (b) Same evolution as (a) but with fixed renormalized charge (e 2 ^ = 50), 
i.e., with the bare charge rescaled with cutoff according to equations (5.9) and (5.10). 
All three evolutions now fall on top of each other. Corresponding to the cutoff values 
of 64, 48, and 32, the values of Z e are .12, .18, and .26, respectively. 



Figure 4: Born Diagram for X exchange which defines the renormalized quartic cou- 
pling Xr in the 1 /N resummed $ 4 field theory. 



The same strategy may be followed at any order of the 1 /N expansion provided 
that the A dependence of the bare couplings are known to the same order in the covari- 
ant treatment. This is an enormous technical simplification over standard adiabatic 
expansion methods for removing the cut-off dependence explicitly. 

In A$ 4 theory the general renormalization algorithm beyond lowest order requires 
knowledge of the (j) 2 X vertex and ^-wavefunction renormalization constants Z\ and 
Z 2 respectively, as well as the coupling renormalization Z\. The first two of these are 
defined in terms of the effective action by the conditions, 

= -i^ki) = i e ' Rl{x ' yh (5 - n) 

with Tr the renormalized proper vertex part, and Q^ 1 the renormalized inverse cj) 
propagator in the case the mean field (j> = 0. In the case (j> ^ we may continue to 
define Q^ 1 as in (5.11) with the understanding that the two-point function for the 
fields (j) and X becomes a matrix with off diagonal (i.e., <j>X) components, so that, in 
particular the renormalized cj) propagator is no longer given simply by inverting Qr 1 ■ 
The other diagonal two-point vertex is 

8 2 S 

V^fx.y) = - ———— = V^(x,y) , (5.12) 

which requires no renormalization, i.e., it is already RG invariant, since it can be 
related directly to the physical Born scattering amplitude as discussed in the third 
paper of Ref. [16] and represented in Fig. 4. Indeed, the renormalized A coupling is 
specified by the Schwinger-Dyson equation for the X inverse propagator, 

AkV) = -V-\q) = Y K -\t^y^ + ( l) r (P.P + S(p) , (5.13) 
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Figure 5: Exact Schwinger-Dyson equation for the X inverse propagator £> _1 . Here 
Q represents the exact unrenormalized (j) field propagator and T the exact unrenor- 
malized vertex function. For the case of QED, the photon inverse propagator has a 
similar graphical representation where the first term is replaced by d~ x and the vertex 

r by iv 

illustrated in Fig. 5. This defines the coupling renormalization Z\(A,m = \/—q 2 ) in 
eq. (5.4) to any order of the 1/N expansion. To lowest order the vertex function, 
T = r = 1 and the result of the previous discussion leading to eq. (5.4) is recovered. 

Now, from the definitions (5.11) it follows immediately that 

= ^-T R (x,y;z) (5.14) 

since the renormalized vertex part Tr may be defined by the X variation of the 
renormalized self-energy appearing in ^ 1 (a3,y). Comparing the last relation with 
the first member of eq. (5.11) implies the Ward Identity, 

Z 1 = Z 2 . (5.15) 

Hence the vertex renormalization Z\ may be computed from the self-energy renor- 
malization Z 2 directly, just as in QED. Indeed, the definitions analogous to (5.11) in 
the QED case are 

oip{x)dip{y) Z/2 

The Ward Identity (5.15) follows immediately for the same reason as before. The 
general polarization tensor is guaranteed to be proportional to q 2 g >il ' — q^q" by the 
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Ward-Takahashi identity following from gauge invariance of the effective action, so 
that 

q,V-\ q r = 0, (5.17) 

and the renormalized charge may be defined by the Schwinger-Dyson equation for the 
photon inverse propagator, 

v-\ q r = -TT-i^V"' - ?v) 

e R\ ( l ) 

= d-\ q y»-zJ A -0^t?{i»g(p)T»(p,p + q )g(p + q )} , (5.18) 

illustrated again by Fig. 5, with the point vertex now representing c? _1 ^(g) and V 
replaced by V 1 . This defines the coupling renormalization Z e (A,m = \/— <? 2 ) to any 
order of the 1/N expansion. To lowest order the vertex function, = Tq = 7 M and 
the result of the previous discussion leading to eq. (5.10) is recovered. 

If there is no mean value of the Dirac field, the Ward Identity guarantees that 
wavefunction and vertex renormalization drops out of the Maxwell equations for the 
mean field and that the charge renormalization Z e (A,m) alone is sufficient to arrive 
at cut-off independent time evolution for the potential. To prove this we make use 
of several facts from the standard covariant analysis. First, we recall that the inverse 
propagator and vertex functions satisfy the Schwinger-Dyson equations, 

g-\q) = G-\q) -^f Q(p + q) P>, p + q) V^p) 

r"(p, q) = T - J ^tr {g(p + r )T fl (p + r, q + r)K:(p + r, q + r,p)g( q + r)} , 

(5.19) 

where K, is the two particle irreducible scattering kernel. These equations are illus- 
trated in Figs. 6. The diagrams contributing to C/ -1 and T to first order in 1/N are 
illustrated in Figs. 7. The advantage of introducing the kernel K, into the discussion 
is due to the fact that gg)C is RG invariant, viz., 

g o k o g = g R o k r o g R (5.20) 

in the condensed notation of Section 2. This means that we may eliminate the bare 
pointlike vertex 7 M in favor of fully renormalized dressed quantities via 

T = + Tr {f* o g o k o g} = i- (rfc + Tr {r£ o g R k r o g R }) . (5.21) 

^1 

Then varying the exact Maxwell equation, 

d v F^ = -ie 2 Ti{Tg[A}} (5.22) 
34 




Figure 6: Exact Schwinger-Dyson equations for (a) Q and (b) T M respectively. 




Figure 7: First two terms in the large N expansion of (a) Q and (b) T M respectively. 
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with respect to the mean potential and using the definitions (5.16) yields 

= iTr{-y»goT v og}8A v 

= zf^Tr {(T R + Tr {T R o^olC fl o Q R }) feo^o Q R } 8A V 

= iTi{{T R + ?i{T R oQ R oK R oQ R })Q R oT R oQ R }6A v , (5.23) 

by eqs. (5.16), (5.15) and (5.21). The point of this exercise is that the Ward Identity 
guarantees the cancellation of wavefunction and vertex renormalization constants 
from the gauge invariant Maxwell equation for the mean potential. Since Z\/Z 2 = 1 
have cancelled out, the remaining loop integration over the skeleton diagrams of the 
last form of eq. (5.23) can have only log divergences which must be absorbed by 
charge renormalization alone. Indeed, comparing the second line of eq. (5.23) with 
the Schwinger-Dyson equation which defines the renormalized charge (5.18) shows 
that this is precisely what happens, since 

V~ lia, 6A v = V R ltlv 6A v (5.24) 

is the RG invariant linear response equation. Setting this quantity to zero allows us 
to study the dynamics of small perturbations in any given background mean field in a 
gauge covariant way independent of the cut-off. Since higher functional derivatives of 
the Maxwell equation involve no divergences whatsoever by simple power counting, 
we have proven that charge renormalization alone is sufficient to render the Maxwell 
equation fully RG invariant, and that the logic of the lowest order renormalization 
carries through to arbitrary order in 1/N without modification. One simply needs 
the corresponding expression for Z e to the given order of 1/iV, which has been given 
in the literature [20]. 

The analogous argument for the the A$ 4 case encounters two complications. First, 
X contains quadratic divergences which must be handled correctly before applying the 
Ward Identity (5.15), and second, the mean field equation for <j> will introduce the need 
for (j) field renormalization over and above the coupling renormalization encountered 
in QED with zero mean Dirac fields. The first complication is handled by shifting 
from X to the renormalized mass, defined by 

Q~\V 2 = -m 2 R ) = . (5.25) 
To first order in 1/N (with </> = 0) we have 

m\ = X + ^S(-r^) = f " ^ f J^G(P 2 ) + ^i-m R ) . (5.26) 
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By definition m 2 R is cut-off independent and RG invariant. Shifting to the physical 
mass pole is necessary in order to remove quadratic divergences on internal lines at 
order 1/N and higher. In the real time formalism the role of m 2 R is played by the 
time dependent effective mass, 

m\t) = X(t) - |- /* dt' f [dk] lm{G(t, t'; k)D(t, t'; k)} (5.27) 
N Jo J 

Since any potential quadratic divergence in m 2 (t) is cancelled by the time indepen- 
dent bare mass counterterm /z 2 , the quantity, m 2 (t) — m 2 (0) is free of such diver- 
gences. At lowest order we were able to show this explicitly by eqs. (5.2)-(5.4). 
At next (and higher orders) the proof is a bit more involved since the quadratic 
divergences appear at the endpoint of the the time integration t' — > t due to the 
singular nature of the integrand(s) in the coincident limit. Hence the divergent lo- 
cal counterterm as t' — > t must be identified explicitly and removed before sensible 
results can be obtained. Since this counterterm is time independent it can be ex- 
tracted most simply by adding and subtracting from the time dependent kernel, 
Xm{G{t,t'; k)D{t,t'; k)} the same quantity evaluated at constant mass, which is the 
same as in the covariant treatment and may be calculated analytically. The differ- 
ence Xm{G{t, t'\ k)D{t, t'\ k) — Go(t, t'\ k)D (t, t'\ k)} then has no quadratic divergence, 
while the remainder has a quadratic divergence structure which is known, and may 
be removed explicitly by integration by parts with respect to t' . The upper limit of 
this integration by parts (at t' = t) gives the local time independent quadratic cut-off 
dependence that is cancelled by subtracting m 2 (0), while the lower limit (at t' = 0) 
gives a finite cut-off independent term which oscillates rapidly (with the cut-off fre- 
quency) for small t. This high frequency "ringing" is a transient result of the "kick" 
to the system coming from our sharp initial condition at definite initial time which 
introduces high frequency components in the Fourier transform, and has been found 
in previous studies of quantum Brownian motion [4]. 

With these prior modifications to eliminate quadratic divergences explicitly, the 
argument leading to eq. (5.23) may be carried over to the scalar $ 4 case line for 
line, with the result that the evolution equation for m 2 (t) — m 2 (0) is completely RG 
invariant, the remaining logarithmic divergences being absorbed by the A coupling 
renormalization alone with Zi/Z 2 = 1 having cancelled from the expression, precisely 
as in the QED case. 

The second difference of the $ 4 interaction from QED arises because of the ex- 
istence of the non- vanishing <j> mean field. Although the renormalization constants 
Zi and Z 2 drop out of the mass equation due to the Ward Identity (5.15), when the 
mean field <j> is non- vanishing then it must be renormalized, so Z 2 {= Zi) will appear 
in the <j> mean field equation. In the covariant formulation Z 2 is the logarithmic wave 
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function renormalization constant given in terms of the derivative of the self-energy 
function on mass shell, namely 



1 dQ~\v 2 



Z 2 dp 2 



2> 



,\ N dp 2 



(5 ' 28) 



Since the field equation for the mean field involves Q~ x o (j) and Z 2 ~ X is a power series 
in 1/AT, the mean field equation is renormalized by multiplying only the first (order 
1/N°) term of the (j) eq. of motion (4.33) by Z 2 ~ X ■ That is, 

+ m 2 (t)) #t)i_ = -l7 fdt'im-m) t Vd*Xm[G>(M'; *)£>(*.*';*)] 
J Zj 2 it N Jo Jo 

(5.29) 

with Z 2 X given by the previous relation is the correctly renormalized equation of 
motion for the (j) mean field to order 1/AT, with the A dependence in Z 2 just can- 
celling that of the momentum integral in (5.29) to first order in 1/N. Naturally, if 
non-zero mean charged fields are considered in electrodynamics they will have to be 
renormalized in the same way. 

One additional remark about the A$ 4 theory and its renormalization is in order. 
It is well known that this theory is trivial, in the sense that the cut-off cannot be 
removed to infinity without vanishing renormalized A. From the point of view of 
practical numerical calculations this is irrelevant since a cut-off will always appear in 
the computer implementation, and the theory with finite cut-off is well defined and 
nontrivial. However, a necessary consequence of this point of view is that the cut-off 
cannot be removed in principle, and one should work in the range where the cut-off 
is not too large, or more precisely where 

< 1 - Z X (A, m) = -^-\ R {m 2 ) In (-) + O (\ R {m 2 ) In (-) *) < 1 . (5.30) 
loir' 1 \mj \ \mj J 

That this is the right condition on the cut-off may be seen either from triviality 
considerations or the senselessness of the theory when the Landau pole is reached 
at Z\ = 1. Similar considerations presumably apply to QED as well, where the 
very weak coupling still affords an enormous range of momenta before the Landau 
pole is reached. So it is always possible to satisfy (5.30) and still have A much 
greater than all frequencies of interest in the non-equilibrium time evolution of the 
fields, if the coupling is small enough. This necessary limitation on the theory has a 
positive side. If the quantity in (5.30) is small it is then permissible to develop the 
expression for D in eq. (4.16) in a power series in A, rather than solving this integral 
equation numerically, thereby recovering its ordinary perturbative expansion. This 
leads to an enormous economy of computer memory since the integral equation for 
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D involves very big arrays and is extremely memory intensive. The evolution of 
the (j) propagator and mean fields is still treated in the full 1/N expansion without 
modification. Further details of this procedure will be presented when we turn to 
numerical methods in future publications. 

Finally, for many applications it is useful to have the energy-momentum tensor 
following from the quantum effective action by variation with respect to the metric 
of spacetime, 

T„ = "^=#^ • (5-31) 

Since general coordinate invariance is maintained in the effective action, the energy- 
momentum tensor is conserved. This is an important property in nonequilibrium 
dynamics which is easily lost if one makes uncontrolled stochastic assumptions or 
approximations in a transport formulation. In the effective action approach, on the 
contrary, conservation of T M „ is automatic, provided the cut-off procedure does not 
violate coordinate invariance. 

In the scalar theory, after scaling out a factor of N we have 

+ £(g^n-dM$ 2 (5.32) 

in terms of the original quantum fields, where the last term proportional to the 
arbitrary parameter £ may be added without affecting conservation. By taking the 
expectation value of this quantity, we may express the conserved energy-momentum 
tensor in terms of the propagator and vertex functions introduced above, viz., 



T^x) 



= T£0,X](a) - i (8 a X - d a d',g(x,x')\ x=x , 

~2N 9 ^ I d?y d?z V(z',x)g(x,x')T(x',y';z')g(y',x) 

where Tj^,[</),X] is (5.32) evaluated on the mean fields, and it must be recalled that Q, 
the full connected two-point function for the $ field is not the inverse of Q~ x defined 
in eq. (5.11) when (j) ^ 0. 

For QED the analogous expression is 

TM = t£[A] - f v (|X tr {y v) g(x, x)} - ^C"°d a d'eV pa (x, x')\ x=x , 

-1 J d 4 x' d 4 y'd 4 z'V a ^(z',x)tr{l v) g(x,x')T a (x',y';z')g(y',x)} 

(5.34) 
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where 

C" 7 = I^C^^" -^^5 + - + *£x<)f - *cL^)^ p ) ( 5 - 35 ) 

and 

Ti[A] = O+rdaAfaK = -Ig^F^Faf, - F^F av (5.36) 

is the stress tensor of the Maxwell mean potential. The next two terms on the first 
line of (5.34) are the contributions to the energy-momentum tensor of the fermions 
(moving in the mean potential A) and photons respectively, while the last term is the 
contribution of the interaction between them which appears first at order 1/N. 

Each energy-momentum tensor contains terms with quartic and quadratic depen- 
dence on the ultraviolet cut-off A. From general coordinate invariance of the effective 
action, these divergent contributions to T M „ must be proportional to the metric of 
spacetime which is flat here. Thus, these divergences may be isolated and re- 
moved rather easily by subtracting from the full T M „ above the same quantity in zero 
mean field(s). The resulting subtracted T M „ is still conserved and now completely 
finite. Indeed, the argument that the current appearing in the Maxwell equation 
(5.22) must be cut-off independent, provided that the bare charge is rescaled with 
the cut-off while keeping the renormalized charge fixed, may be taken over to the 
subtracted energy momentum tensor as well. Like the mass in the $ 4 theory, once 
the power law divergences are removed from T M „ the resulting quantity is RG invari- 
ant, and may be interpreted as the physical energy-momentum of the nonequilibrium 
field theory evolution. Since this T M „ is computed from the same propagator and 
vertex functions appearing in the evolution equations, the pressure, energy density 
and transport characteristics of the QED plasma or quantum $ 4 theory may be stud- 
ied, and useful information about the approach to hydrodynamic behavior and/or 
an effective equation of state obtained. The effective equation of state for the QED 
plasma to lowest order in 1/N has been discussed in Ref. [10]. Detailed numerical 
results for these quantities and further applications of the large N CTP method are 
in preparation and will be presented in future publications. 

6 Conclusions 

In this paper we have presented a general approach to the nonequilibrium evolution 
of a closed system of quantum fields. The large N expansion permits a clean separa- 
tion into mean fields and their fluctuations and constitutes a controlled approxima- 
tion scheme to the infinite tower of coupled Schwinger-Dyson equations of quantum 
field theory. We have derived the equations from an effective action principle which 

40 



preserves all classical symmetries (which are not anomalous) at the quantum level. 
This is important because of the central role invariances play in the renormalization 
procedure through the Ward Identities, which are in danger of being obscured in a 
noncovariant time evolution problem. Nevertheless, manifest causality of the time 
evolution is enforced by the Schwinger-Keldysh CTP formulation of the effective ac- 
tion principle. We have sketched a renormalization procedure involving an ultraviolet 
cut-off which is well-suited to numerical solution of the equations on a computer, and 
demonstrated explicitly the practicality of this scheme at lowest order of the 1/N 
expansion. 

For definiteness, throughout the paper we have developed the general approach 
in the framework of two particular and familiar quantum field theories, viz., A$ 4 
and QED. From this beginning there are three well-defined vectors for future work. 
The first consists in carrying out practical numerical computations in these theories 
for specific applications, for example, to the evolution of disoriented chiral conden- 
sates in heavy-ion collisions or to e + e~ particle production and shorting of strong 
fields in astrophysical plasmas. These applications are clearly interesting in their own 
right. The second direction to pursue is the use of these particular field theories as 
model systems for the study of more general phenomena, such as dissipation and 
decoherence in closed quantum systems. The emergence of an effective Boltzmann 
or transport equation description from fundamental time reversal invariant dynamics 
may be studied in these realistic field theories in a controlled way without additional 
stochastic assumptions. Finally, because of the existence of a gauge invariant action 
principle, the general method followed in this paper can be extended to non-abelian 
gauge theories such as QCD and gravitation without essential difficulty. This will 
make it possible to take into account consistently the backreaction of quantum fluc- 
tuations on the nonequilibrium evolution of a mean color or metric field and open up 
interesting applications in studies of the quark-gluon plasma, black hole decay, and 
cosmological models. We plan to take up each of these lines of research in subsequent 
publications. 
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